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Abstract 



Preliminary results are presented on the gluon density at low x obtained from a QCD 
I analysis of ZEUS 1994 F2 structure function data combined with those from NMC. Also 

given are estimates of the experimental error on the e^p NC Born cross section at large 
X and Q^. This estimate is obtained from propagation of the statistical and systematic 
Q^ ! errors on fixed target structure functions. 
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P ■ Introduction 

On , 

In perturbative QCD the scaling violations of the F2 structure functions are caused by gluon 
bremsstrahlung from quarks and quark pair creation from gluons. In the low x domain accessible 
gn! at HERA the latter process dominates the scaling violations. In this report I present preliminary 
r-| ! results on the gluon momentum density extracted from a NLO QCD anlysis of F2 structure 
functions measured by ZEUS in 1994. 

^ ■ With the present integrated luminosity of about 20 pb^^ at HERA the region of large x (~ 0.5) 
^ . and (~ 10^ GeV^) becomes accessible. Both ZEUS lU and HI have observed an excess 
of events compared to the standard model predictions in this hitherto unexplored region. The 
uncertainty in these predictions is dominated by that on the parton distributions in the proton. 
I will summarise the results of a QCD analysis by ZEUS of fixed target structure function data 
at high X (> 0.1) yielding an estimate of the experimental error on the NC Born cross section 
of e~^p deep inelastic scattering at HERA. 



QCD analysis of ZEUS 1994 data 

The data used in the fit were the ZEUS 1994 nominal vertex data which cover a kinematic 
range of 6.3x10"^ < x < 0.5 and 3.5 < < 5000 GeV^ together with the low measurements 
(shifted vertex, ISR) [Q. The latter datasets extend the kinematic range down to x = 3.5 x 10~^ 
and = 1.5 GeV^ albeit with larger statistical and systematic errors. NMC data on and 
F2 constrain the fit at high x. To remove possible contributions from higher twist effects at 
large x the NMC data below = A GeV^ were discarded. 

^Presented on behalf of the ZEUS collaboration at the 5th International Workshop on Deep Inelastic Scat- 
tering and QCD, Chicago, IL, April 14-18, 1997 
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At the input scale Ql = 7 GeV^ the gluon distribution ixg), the singlet quark distribution (xS) 
and the difference of up and down quarks in the proton {xAud) were parametrised as 

xg{x,Ql) = AgX^^{l-xy^{l+^gx) 

xJ:ix,Ql) = A,x^^{l-x)^^{l + esV^ + 7sx) (1) 
xAUx,Ql) = A^,x'"'{l-xy- 

The strange quark distribution was assumed to be 20% of the sea at = 4 GeV^ 0. The sea 
quark density was otained by subtracting the valence (taken from MRSD'_) from the singlet 
distribution. The gluon normalisation, Ag, was fixed by the momentum sumrule. The input 
value for the strong coupling constant was set to the result of ref. [0: as(M|) = 0.113. 

The input parton distributions were evolved in NLO in the MS scheme with / = 3 light flavours. 
The charm contribution to the F2 structure function was calculated in NLO from the evolved 
distributions as described in |^ with the charm mass set to 1.5 GeV. Contributions from bottom 
are estimated to be small and were neglected. 

In addition to the 11 parameters describing the parton distributions one (two) normalisation 
parameters for the ZEUS (NMC) data were left free in the fit. In the computation of the 
only statistical errors were taken into account. For each dataset the quantity [(A^ — 1)/AA^]^ 
was added where AA^ is the quoted normalisation error of the dataset. 

The fit yielded a good description of the data as shown in Fig. |I[ Adding the statistical and 
systematic errors in quadrature the = 463 for 408 datapoints and 13 free parameters. The 
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Figure 1: The structure function versus for fixed values of x. The errors shown are 
the statistical and systematic errors added in quadradure. The solid curves correspond to the 
NLO QCD fit described in the text. The dashed curves show F2 without the contribution from 
charm. 
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fitted normalisation parameters were 97.4% for ZEUS and (99.3%,98.9%) for the NMC 90 GeV 
and 280 GeV datasets respectively which is well within the quoted normalisation errors. 

Fig. ^ shows the gluon momentum denstity obtained from the fit at the input scale Ql = 1 
GeV^. The inner (outer) shaded band indicate the statistical and the statistical and systematic 
error added in quadrature. It is seen that at the lowest value of x = 3 x 10~^ the gluon is 
determined with an accuracy of about 20%. At x = 4 x 10~^ the total error is ~12% which is 
a large improvement compared to the ~50% error on the gluon distribution obtained from the 
QCD analysis of the ZEUS 1993 F2 data At the input scale the momentum fraction carried 
by quarks (gluons) was found to be 0.555 (0.445). The fitted values of the parameters describing 
the low X behaviour of the quark singlet and the gluon distributions are 5s = — 0.23±0.01±0.02 
and 6g = —0.24 ± 0.02 ± 0.05 where the first error is statistical and the second systematic. 

The gluon distribution evolved to = 20 GeV^ is shown as the shaded (total) error band 
in Fig. 0b. Also plotted are the gluon densities from the recent parton distribution sets 
MRS(Rl) 0, CTEQ4M Jill and from GRV(HO) p[. Whereas the agreement with MRS 
and CTEQ at low x is excellent, the steep gluon obtained from the dynamical evolution by 
GRV is inconsistent with the result of this analysis. 




Figure 2: The gluon distribution as function of x for a fixed value of = 7 GeV^. The 
inner (outer) shaded band corresponds to the statistical (total) error. Also shown is the gluon 
distribution obtained from the ZEUS 1993 analysis, (b) The gluon distribution from this 
analysis (shaded band) at = 20 GeV^ compared to those from MRS(Rl), CTEQ4M and 
GRV(HO). 
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Extrapolation of high x structure functions to large 



In this section we describe a QCD analysis of high x structure function data. The aim of this 
analysis is to estimate the experimental uncertainty on the ep neutral current (NC) cross-section 
at large x ~ 0.5 and ~ 3 x 10^. This error estimate is used in [|l| to judge the significance 
of a possible excess of events observed at HERA in this region [|T], 0. 

The data used in the fit were the proton and deuteron F2 data from SLAC [113, BCDMS [IH 
and NMC |15| together with xF^^ from CCFR Since we are only interested in the high 



X domain, data below x = 0.1 were discarded. Also applied were the cuts > A GeV^ and 
> 10 GeV^ to remove target mass and higher twist effects. F2 data above x = 0.7 were 
discarded to eliminate possible contributions from Fermi motion in deuterium. 

The QCD analysis of these data is similar to that presented in the previous section except (i) the 
charm (bottom) quark distribution was generated dynamically from the threshold Qc(fe) ~ ^ 
(25) GeV^ (ii) no momentum sumrule was imposed (iii) all normalisations were kept fixed to 
unity and (iv) the valence quark distribution, constrained by the CCFR xF^ data, was left free 
in the fit with a functional form identical to that of xA^d in Eq. (|l]). From the evolved quark 
distributions the F2, Fl and xF^ structure functions in ep scattering were calculated in NLO 



including contributions from Zq exchange and jZ interference 



The experimental systematic errors were propagated to the covariance matrix of the fitted 



parameters using the technique described in |18|. In total 24 independent sources of systematic 



error were included taking properly into account the correlations between the systematic errors 
of the NMC datasets. 

It is convenient to express the NC cross section for e^p scattering as 

where y = Q'^/xs, s ~ AEeEp is the square of the ep centre of mass energy and Y± = 1±(1 — y)^. 

In Fig. ^ we show the cross section for e~^p scattering at HERA energies {E^ = 27.5 Ep = 820 
GeV) as a function of for fixed values of x. The strong decrease in the cross section above 
Q"^ ~ 5000 GeV^ is due to the contribution from xF^. At high x ~ 0.5 and ~ 3 x 10^ GeV^ 
the error on is estimated to be about 9%, see Table . The estimated errors include (small) 
contributions from an assumed uncertainty of Aas(M|) = 0.005 and an error of 50% on the 
strange quark content of the proton. 

In Fig. the cross sections obtained from this analysis are compared to those calculated from 
the parton distribution sets MRSA and CTEQ3M |2^. It is seen that the differences 
are much smaller than the experimental errors: comparison of results obtained from different 
parton distribution sets do not yield a reliable estimate of the uncertainty in the cross sections. 

(GeV^) X 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
1000 4444458 15 31 

30000 11 9 8 9 14 30 



Table 1: The relative uncertainty in a'^^Q in percent. 
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Figure 3: (a) The cross-section a (see text) for e+p scattering at Hera energies, (b) The result 
of this analysis compared to cross sections calculated with the parton distribution sets MRSA 
and CTEQ3M. 
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